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Introduction.  -We  consider- ^stabilization  (minimum  time  to  the  zero  solution)  of  wave  and 
beam  equations  via  endpoint  velocity  control.  In  general,/  consider  a  control  system  described 
rv  a  partial  differential  equation,  with  dependent  variable  denoted  w  =  w(t,x),  which  admits 


V  =  — —  as  a  symmetry  generator  and  hence  a  conservation  law  which  implies  conservation  of 

£7t 


energy,  denoted  e(w(t,*)).  For  the  problems  considered,  it  is  shown  that  — —  e(w(t,'))  is  a  con- 

dt 


vex  function  of  the  control  u;  the  value  u*  which  minimizes  this  can  be  expressed  in  feedback 
form,  i.e.,  if  w(t,-)  belongs  to  a  Banach  space  B  and  /  :  B  — *  is  a  closed  linear  map  (a 
finite  dimensional  measurement  on  the  state  of  thi  system;  the  output  of  sensors),  then  the 
optimal  control  will  have  the  form  u*  =  f(w(t,-)). 

The  energy  decay  of  solutions  of  wave,  or  beam,  equations  via  boundary  control  has 
received  much  attention  in  the  recent  literature.  In  111,  the  wave  equation 


/  %  d2\v  ^  d  ,  ,  dw  n 

P(x)  TT"  ~  £  K(*))  T  =  0  > 


<9t2  dx,  v  "  5Xj 


(1) 


with  t  >  0,  x  €  C  lRn,  n  open,  bounded,  connected  and  with  boundary  dQ  =  T,  is  studied. 
Assume  p(x)  >  p0  >  0;  the  matrix  A  =  (Qr(J(x))  is  positive  definite  symmetric  (i.e.,  the  inner 
product  (£,  A(x)£)  >  o0j  for  some  o0  >  0).  Partition  the  boundary  as  T  =  T,  U  with 
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rt  relatively  open  and  require  w(t,x)  =  0  on  w(t,x)  =  h(t,x)  a  control  for  x  £  [\>.  The 
energy  of  a  solution  w  of  (l)  is 

e(w(t,  •))  =  (-J-)  /  (p(x)w2(t,x)  +  (wx(t,x),  A(x)wx(t,x))]dx  .  (o) 

-  n  w 

If  w  satisfies  (1)  and  t^x)  denotes  an  outward  normal  to  dCl,  the  divergence  theorem  yields 

~  e(w(t,  •))  =  /  wt(t,x)(wx(t,x),  A(x)^(x))dx  .  (3) 

Change  control  from  "position  h"  to  velocity  u(t,x)  =  ht(t,x)  =  wt(t,x),  x  6  rt.  From  (3),  a 
natural  choice  to  decrease  energy  is 

u(t,x)  =  wt(t,x)  =  ~ (wx(t.x),  A(x)i^x)),  xCTj.  (4) 

A  great  deal  of  the  literature  is  devoted  to  an  analysis  of  (3)  with  control  choice  (4).  In  [I]  an 
invariance  principle  (extended  LaSalle  principle  as  given  by  Hale)  is  used  to  show  e(w(t,  •))  —  0 
as  t  — *•  00.  Under  more  special  conditions  on  the  shape  of  f2  and  with  p{\)  —  1,  A(x)  =  I,  the 

convergence  rate  e(w(t,-))  =  c  t  >  0,  is  obtained,  where  c  depends  on  the  initial  condi- 

tions.  It  is  well  known  [2j,  [3,  Thm.  1],  that  solutions  of  the  one-dimensional  wave  equation  can 
be  driven  to  zero  in  finite  time,  conditions  which  insure  this  for  quasi-linear  hyperbolic  equa¬ 
tions  can  be  found  in  [4],  For  more  general  systems,  exponential  asymptotic  estimates  of  rate  of 
energy  decay  can  be  found  in  [5],  [6].  Recent  results  involving  the  physical  character  of  sensors 
and  actuators  in  energy  decay  are  given  in  [7j.  Energy  methods  for  boundary  control  of  the 
beam  equation  can  be  found  in  (8),  (9j. 

Our  purpose,  here,  will  be  to  stress  geometry,  optimality,  the  feedback  nature  of  the  con¬ 
trols  (with  it  the  sensor  placement  problem)  and  to  suggest  interesting  areas  for  further 
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1.  The  Endpoint  Control  of  a  Vibrating  String 

Let  w(t,x)  denote  the  vertical  displacement  at  position  x,  0  <  x  <  1,  and  time  t  of  a 
string  of  length  1.  The  equation  of  motion  is  taken  as 

wu  =  wxx-  0  <  x  <  1,  t  >  0  .  (5) 

Initial  Conditions  (IC):  w(0,x)  =  wt(0,x)  =  r?(x)  (6a) 

Boundary  Conditions(BC):  w(t,0)  =  0,  wt(t,l)  =  u(t)  =  control  (6b) 

The  problem  is  to  use  the  velocity  control  (or  one  could  use  height  or  right  endpoint  position 
control)  to  drive  the  solution  to  w^.x)  =  wt(t1(x)  =0,  0  <  x  <  1,  for  minimum  time  tj.  It  is 
well  known,  [2],  [3],  that  for  the  length  one  string  and  arbitrary  initial  data  this  can  be  achieved 
in  time  tj  =  2  and  no  smaller  time  suffices.  Indeed,  calculations  give  the  explicit  formula 

'-(jJW-O-iK!-*)).  o  <  t  <  i 

u*(t)  =  '  (7) 

-(yX^t-i)  +  *(t-i))f  i  <  t  <  2 

For  this  simple  problem,  we  see  in  (2)  that  Tj  =  {0},  T2  =  {l},  while  if  we  denote  the 
solution  of  (5),  (6)  as  w(t,x,u)  to  explicitly  show  the  u  dependence, 

~  e(w(t,-,u))  =  u(t)wx(t,I,u)  .  (8) 

A  natural  choice  (as  in  (4))  is  to  choose  u(t)  =  — kwx(t,l,u),  k  >  0.  The  interesting  thing  (as 
will  be  shown)  is  that  the  right  side  of  (8)  is  quadratic  in  u  and  has  minimum  at  k  =  I,  i.e., 
u*(t)  =  — wx(t,l)  is  the  optimal  control.  Indeed,  one  can  see  this  by  considering  (5)  with  initial 
conditions  (6a)  and  (nonstandard)  boundary  conditions  w(t,0)  =  0,  wt(t,l)  =  — wx(t,l).  This 
problem  can  be  explicitly  solved  by  standard  techniques;  one  readily  fin^s  that  the  solution 
satisfies  wt(t,l)  =  u*(t)  as  given  by  (7). 
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Finally,  for  the  endpoint  controlled  string,  we  have 

Proposition  1.  The  control  choice  u*(t)  which  pointwise  minimizes 

-7-  e(w(t,  •,  u))  =  u(t)wx(t,  1,  u)  satisfies  u*(t)  =  — wx(t,  1,  u*)  and  is  the  optimal  control, 
dt 

To  prove  this,  solve  (5)  with  initial  conditions  (6a)  and  boundary  conditions  (6b),  leaving 
u(t)  arbitrary,  in  a  region  R  containing  the  right  boundary  x  =  1.  Explicitly,  for 
R  =  {(t.x):  x  +  t>  1,  x  —  t>0,  x  <  1}  one  may  show 

2-x-l  x+t-i 

w(t,x,u)  =  “{~-)<^(2— x— t)  4-  (~)^(x—t)  4-  (^-)  /  n(s)ds  +  <5(1)  4-  /  u(s)ds  . 

This  allows  us  to  compute  lim  wx(t,x,u),  which,  for  0  <  t  <  1,  is 

x—i- 

wx(t,l,u)  =  <^(1— t)  —  t](  1— t)  4-  u(t).  Substituting  this  in  the  equation  for  the  time  derivative  of 
e(w(t,*,u))  shows  e(w(t,-,u))  -  u(t)[tf'(l-t)  -  r/(l-t)  4-  u(t)],  1  <  t  <  1.  Let 

f(u)  =  u(<y— t/+u).  This  quadratic  has  a  minimum  at  u*(t)  =  — (~)(d>'(l—  t)  —  //(l  —  t)), 

0  <  t  <  1.  A  similar  analysis  shows  that  for  1  <  t  <  2,  u*(t)  =  — ( J-)(<^,(t — l )  4-  r/(t— 1)). 

2 

Comparison  with  (7)  shows  u*  is  optimal. 

The  point  to  be  stressed  is  the  quadratic  nature  of  the  right  side  eq.  (S)  which  assures  an 
absolute,  pointwise,  minimum.  It  is  also,  geometrically,  interesting  to  note  that  the  intuitive 
idea  of  choosing  the  right  endpoint  (i.e.,  wt(t,l))  to  be  the  negative  of  the  slope  of  the  wav.'  at 
this  endpoint  (i.e.,  — wx(t,l)),  is  optimal. 

2.  Endpoint  Energy  Control  of  a  Beam 

For  ease  of  exposition  we  use  the  overly  simple  beam  model 


p wtt  +  EIwxxxx  =  0,  0  <  X  <  /,  t  >  0  (g; 

I.C.:  w(0,x)  =  <^x),  wt(0,x)  =  r?(x),  0  <  x  <  l 

We  assume  a  free  right  end  yielding  right  boundary  conditions 

BC(R):  wJt.O  =  wxxx(t,0  =  0 

The  left  boundary  conditions  will  involve  the  controls  and  yield  three  natural  problems 
Before  discussing  these,  we  note  the  energy  now  is 


eNt,-))=(|)/(pwt2  +  EIw^)dx 
z  o 


which,  using  B.C.(R),  satisfies 


—  e(w(t,  •))  =  Elwl(t,0)wxxx(t,0)  -  EIwxl(t,0)wxx(t,0)  .  (H) 

Problem  1.  The  left  endpoint  velocity,  i.e.  wt(t,0),  is  the  control  and  this  end  of  the  beam  is 
held  so  as  to  always  have  wx(t,0)  =  0.  This  eives  left  boundary  conditions 

B  C  (L):  wt(t,0)  =  u(t)  =  control,  wxl(t,0)  =  0 

Problem  2.  Consider  the  left  end  of  the  beam  fixed  at  zero  with  the  initi  tl  slope,  wx(t,0),  the 
control.  More  precisely,  we  take  the  time  rate  of  change  of  the  initial  slope  as  control. 

B.C.(L):  w(t,0)=0  (or  wt(t,0)  =  0),  wxl(t,0)  =  u(t)  =  control 
Problem  3.  Both  initial  (left  end)  height  and  slope  are  controllable. 

B.C.(L):  wt(t,0)  =  u,(t),  wxt(t,0)  =  u.(t). 

Motivated  by  the  results  of  the  one-dimensional  wave  equation,  we  see  each  problem  sug¬ 
gests  a  natural  control  choice  to  decrease  energy,  i.e.,  in  problem  1  choose  u(t)  =  — wxxx(t,0);  in 
problem  2  choose  u(t)  =  Wxx(t,0),  in  problem  3  choose  u,(t)  =  — wxxx(t,0),  u2(t)  =  wvx(t.O).  Are 


6 


these  choices  (or  positive  multiples  of  them)  optimal?  For  the  sake  of  brevity,  we  examine  only 
problem  2  for  which  (11)  becomes 

-j-  e(w(t,-,u))  = -Elwxl(t,0)w  (t,0,u)  (12) 

at 

=  — Elu(t)wxx(t,0,u) 

where  (again)  stress  of  dependence  of  w  on  u  is  made  clear  where  needed.  The  choice 
u(t)  =  w^t.O)  physically  means  pick  the  velocity  of  the  initial  slope  to  be  the  initial  curvature; 
a  reasonable  scheme!  Also,  physically,  one  expects  |u|  sufficiently  large  implies 

—  >'  (t,,,u))>0  while  u  =  0  implies  ■—  e(w(t,  *,  u))  =  0,  i.e . ,  — -  e(w(t,  •,  u))  can  be 

dt  dt 

\pected  to  be  convex  in  u  with  (when  w^t.O.u)  ^  0)  a  negative  absolute  minimum.  To  com¬ 
pute,  explicitly,  one  must  calculate  w^tjO.u).  This  may  be  accomplished  via  the  Heaviside 

d 

operational  calculus,  see  [10,  pg.  51].  Briefly,  let  p  =  — — ,  h(t)  =  .vx(t,0)  so  ph  =  u,  wH 

(J  l 

denote  the  solution  of  (9)  with  homogeneous  boundary  conditions  and  (symbolically)  express  the 
solution  of  problem  2  as  w(t,x,u)  =  wH(t,x)  +  T(x)h(t).  The  nonhomogeneous  part,  i.e.  v  =  Th, 
must  satisfy  (EIT^^  +  /)p2T)v  =  0,  Tx(0)  =  1,  T(0)  =  0,  T^/)  =  0,  T^/)  =  0.  Considering  p 
as  a  parameter,  and  letting  a2  =  pp2/EI,  one  finds 

T(x)  =  Cjex N//^2cos(x  ^a/2  4-  c2ex N//^2sin(x  ^a/2)  4-  c3e_x  V'^2cos(x  ^a/2)  -f  c4e-x  V'n^'sin(x  ^a/2 

We  are  interested  in  Txx(0)h,  from  which  it  follows  only  c2,  c4  are  of  importance.  Using  the 
boundary  data  and  T  as  above,  with  tedious  calculations  one  finds  e2  ~  'v/p,  c.,  ~  'v^p.  Not¬ 
ing  that  or  is  a  constant  multiple  of  p,  Txx(0)h  ~  p3^2h  =  p1/,2u  which  (see  [10!)  implies 


TjO.hft)  ~  y-  /  ds  . 

dt  0  Vt_^ 


Finally,  for  a  computable  constant  r  e  have  the  behavior 
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e(w(t,  *,  u))  ~  u(t)w^(t,0)  4-  cu(t)  — -  /  -  ^i===-  ds  .  (13) 

at  dt  o  V t— g 

Questions:  1.  Can  one  carry  out  the  above  computations  precisely  and  by  choosing  u  to 
minimize  the  right  side  of  (13),  explicitly  calculate  an  optimal  (for  minimizing  time  to  zero  of 
e(w(t,  *,  u))  control  u*? 

2.  Can  one  use  e(w(t,-))  =  y l(t),  w^t.O)  =>'2(0  as  two  of  a  finite  number  (say  n)  of 
measurements  on  a  solution  w  such  that  ylt  •  -  ,yn  satisfy  a  controlled  system  of  ordinary 
differential  equations?  If  so,  one  would  use  control  to  drive  yi  to  zero. 
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